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Abstract 

The competition in the pinning of a directed polymer by a columnar pin 
and a background of random point impurities is investigated systematically 



5t ^ using the renormalization group method. With the aid of the mapping to 

the noisy-Burgers' equation and the use of the mode-coupling method, the 
directed polymer is shown to be marginally localized to an arbitrary weak 
columnar pin in 1 + 1 dimensions. This weak localization effect is attributed to 
the existence of large scale, nearly degenerate optimal paths of the randomly 
pinned directed polymer. The critical behavior of the depinning transition 
above 1 + 1 dimensions is obtained via an e-expansion. 
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I. INTRODUCTION 



The statistical mechanics of an elastic manifold embedded in a medium of random point 



defects has been the subject of many studies in the past decade Such systems are 

encountered in a variety of contexts, ranging from the fluctuations of domain walls in random 
magnets , to the pinning of magnetic flux lines in dirty superconductors |]^-^ . Over 

the years, many theoretical methods have been developed to understand these systems P|-pT|. 
In particular, through a mapping to stochastic hydrodynamics [^^-|T^, we now know many 
properties of one-dimensional manifolds (directed polymers) in random media. There have 



also been numerous numerical simulations; some recent studies can be found in Ref. 18-20 



Recently, strong flux pinning effects |21-23] exhibited by samples of high temperature 



superconductor with extended defects such as columnar faults and twin planes lead nat- 



urally to the investigation of competition between extended and point defects |24,25|. It 



has been argued in the case of many interacting flux lines that pinning by extended de- 



fects are weakened by point defects . There have also been a number of studies on the 
competing effects between extended and point defects on a single flux line or a directed 
polymer pB|-|55|. An early study of this type was done by Kardar who investigated nu- 
merically the pinning of a directed polymer to a line defect in the presence of a background 
of point defects in 1+1 dimensions. The polymer was found to depin from the line defect, 
if the pinning potential of the line defect is smaller than a certain threshold value. Critical 
behavior associated with the depinning transition was later investigated in more detail by 
Zapatocky and Halpin-Healy |2^. The results of Refs. |2^ and |^ have been challenged 



by Tang and Lyuksyutov |28|, who argued against the existence of a depinning transition 
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in two dimensions, based on large scale simulations and approximate renormalizat ion-group 
analysis on a hierarchical lattice. Tang and Lyuksyutov proposed instead that a directed 
polymer is always localized, although weakly, to the line defect in 1+1 dimensions, and that 
the depinning transition only exists above 1+1 dimensions. This conclusion is supported by 
a recent study of Balents and Kardar [^, who performed numerical simulations, and also 
developed a functional renormalization group analysis by describing the directed polymer as 
a generalized D-dimensional manifold and studying the depinning transition using a double 
expansion (in the dimensionality of the manifold D = 4 — 6 and the dimensionality of the 
embedding space d = dc{D) + e.) On the other hand, analysis of the directed polymer itself 
by Kolomeisky and Straley have led to different conclusions when different renormalization 
group Ansatz were used |3^j31|j3^ . 



In this paper, we give a systematic analysis of the competition in pinning between point 
and line defects for a directed polymer. By exploiting known knowledge of the randomly 
pinned directed polymer in 1+1 dimensions in the absence of any extended defects WTl, we 



construct a renormalization-group analysis directly in 1+1 dimensions, the critical dimension 



for this problem. Our results prove the conclusion of Ref. p8|, that the polymer is always 
pinned at and below 1+1 dimensions, and that the depinning transition only occurs above 
two dimensions. The existence of weak localization at the critical dimension is understood 
in terms of the anomolous large scale excitations of the directed polymer, known as the 



"droplet excitations" [|T^. Critical scaling behaviors at the depinning transition are then 
obtained using a 1 + 1 + e dimensional expansion. 

The paper is organized as follows: We first introduce the known properties of the ran- 
domly pinned directed polymer in Section II. We consider the effect of pinning by a line 



defect using phenomenological scaling arguments, which estabhshes 1+1 dimensions to be 
the critical dimension. In Section III, we attempt to solve the problem using the replica 
method at the critical dimension. Despite the existence of an exact solution via the Bethe 
Ansatz [|n],|T^ when the extended defect is absent, we failed to develop a systematic and con- 
trolled way of incorporating a weak extended defect. We next use an uncontrolled Hartree 
approximation, which contains the right physical ingredients, but overestimates the effect of 
point disorder. A systematic analysis taking advantage of known mapping of the directed 



polymer to the noisy-Burgers' equation is presented in Section IV. Through explicit 
calculations, facilitated by the use of the mode coupling approximation in 1 + 1 dimensions, 
we construct the renormalization-group recursion relation and derive various scaling behav- 
iors near the critical dimension. We show that the results can be understood naturally in 
terms of the anomolous droplet excitations [1^. In Section V, we generalize our method to 



describe the directed polymer pinned by a variety of extended and point defects, including 
columnar pin of extended range and trajectory, and point defects with long range corre- 
lations. Some useful relations about the randomly pinned directed polymer, in particular, 
the mode coupling approximation, is summarized in Appendix A, and a number of detail 
calculations are relegated to Appendix B. 
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II. PHENOMENOLOGICAL SCALING ANALYSIS 



We describe the statistical mechanics of a randomly pinned directed polymer of length t 
by the Hamiltonian 
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where ^{z) G 3?'^-^ denotes the transverse displacement of the polymer in d = dj_ + 1 em- 
bedding dimensions (with d± being the co-dimension), k, is the polymer line tension, and 
ri{-x, z) describes a background medium of uncorrelated point defects. The random potential 
is assumed to be Gaussian distributed with mean zero and the variance 



r7(x, z)r]{^, z') = 2Ar^ (x - x.')5{z - z) , (2.2) 

where the overbar denotes disorder average. 

The randomly pinned directed polymer has been the topic of detailed investigations 
during the last decade [Q- [^. The emerging qualitative picture for a long polymer {t — > oo) 
is the following: At low temperatures, the static properties of the polymer are dominated 
by the random potential and are controlled by one or a few optimal path(s) which minimize 



the total free energy [|1J,|T^]. To take advantage of fluctuations in the random potential ?7, 
the optimal path executes large transverse wandering [solid line of Fig. 1(a)]. If we fix one 
end of the polymer, then the root-mean-square displacement of the other end of an optimal 
path is 

X(t) = Bt^, (2.3) 

where 3/4>(^>l/2isan universal exponent which depends on d. The result C, = 2/3 in 
d = 1 + 1 has been obtained exactly by a number of methods [|ll|,|l^. Numerical calculations 



give (^5/8ind — 2 + 1. Displacement of the fixed end of the polymer by a distance 
r <^ X{t) typically result in a rearrangement of the optimal path within a segment of length 
T ~ (r/S)^/^ from the fixed end, as shown in Fig. 1(a). For sufficiently large displacement, 
i.e., for r ^ Bt'^, the optimal path is completely changed [dotted line of Fig. 1(a)], with a 
completely different free energy [Fig. 1(b)]. The typical free energy difference AFq between 
such independent paths is 

AFo ^ At^, (2.4) 

with an identity 

0^2C-1 (2.5) 

relating the two exponents. AFq also sets the scale of the free energy difference between the 
optimal path and a typical path. 

Thermal fiuctuations do not change this picture qualitatively. They merely wipe out the 
fine structure of the random potential, leading to temperature dependent amplitudes. For 
instance, B oc (A//tT)^/^ and A oc kB'^ ind=l + l. Ind<2 + 1 dimensions, this picture 
is in fact correct for all temperatures. In d > 2 + 1, the polymer undergoes a continuous 
phase transition and becomes dominated by thermal fluctuations rather than disorder at 
sufficiently high temperatures. We shall however be focused on the more interesting low 
temperature phase throughout this paper. 

Let us consider the influence of an additional pinning potential, Up{'x), in the form of 
a fine defect with a short range a, located at the origin. We model the defect by the 
Hamiltonian 

n,= fdzUpi^iz)], (2.6) 
Jo 
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with Up{x) — —U for |x| < a, and C/p(x) = for |x| > a. Clearly, such an attractive 
potential favors configurations of the polymer close to the origin, i.e., it attempts to localize 
the polymer. On the other hand, the polymer will try to wander away from the defect in order 
to take advantage of fluctuations in the background random potential r){x.,z). Depending 
on the outcome of this competition, the pinning potential Up may have a number of possible 
effects: It may always or never localize the polymer for any finite pinning strengths U, or it 
may localize the polymer only above a certain critical pinning strength Uc, thus giving rise 
to a depinning transition. 

To gain an understanding of the competing effects, we first consider the pure problem 
with 77 (x, z) = 0. Let 

iy(°)(x,t) = /^"'%[|]e-<'+^^)/^ (2.7) 

^(0,0) 

be the Boltzmann weight of propagating the polymer from (0, 0) to (x, t) , with — 
'Ho[rj — 0] and the superscript (0) indicating the absence of the random potential, then l^*^*^) 
satisfies the diffusion equation 

T|-W^(o)(x,i) = — V'W^(°) - C/p(x)iy(°)(x,i). (2.8) 
(ji 2/^ 

Our problem is equivalent to that of a d^-dimensional (imaginary-time) quantum mechanical 
particle with a pinning well Up at the origin. In the thermodynamic limit i — > 00, the 
statistical mechanics of the polymer is given by the ground state of the quantum problem 
whose solution is well known: There is always a bound state for arbitrarily weak pinning 
potential in d± < 2, where the polymer is always bounded to the line defect. For d± > 2, a 
critical pinning strength of the order Uq*^ = T^/ (2«;a^) is necessary to have a bound state. 
A depinning transition occurs a,t U — U^^ ~ Uq'\ Fig. 2 summarizes the phase diagram of 
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the pure system in various dimensions. 

It is useful to understand these resuhs by simple physical considerations: Assume first 
that the polymer is completely localized within the pinning well. Then the energy gained 
(compared to a free polymer) is U per length, and the entropy cost (or the kinetic energy 
cost of localizing a quantum particle) is of the order Uq'^ per length. Thus ii U ^ Uq'\ 
the total free energy of the polymer is F^^^ — —U + Uq'^ < 0, and the polymer is locahzed. 
However, ii U ^ Uq^^ , then > and the polymer will not be completely localized to the 
well since the free energy exceeds that of a free polymer in this case. An alternative scenario 
when U ^ Uq^^ is to have the polymer completely delocalized, ignoring the existence of the 
weak pinning potential. This is self-consistent for d± > 2, since for a polymer with one end 
fixed at the origin, the accumulated probability of a long free polymer returning to the origin 
is t(2-(ix)/2 _^ g 'jj^^g ggg lY^Q^ ^Qj. ^ 2, the polymer is essentially free for U <C ^7o°\ 
and is completely locahzed if C/ » C/g'^^- The simplest phase diagram is then to have a 
depinning transition separating the pinned and free phases, occuring at C/ = U^^ ~ Uq\ 
On the other hand, for dj. < 2, the probability of a long free polymer returning to the origin 
is of order 1. Thus the completely delocahzed phase is not self consistent here, and the effect 
of a weak pinning potential can never be ignored in the limit t — > oo. 

The above can be stated more quantitatively by comparing the average pinning energy 
{Hi)^^^ with the available thermal energy T. A useful parameter to focus on is the dimen- 
sionless ratio g^^^ = {Tii)'^^^ /T . The polymer is free if g'^'^'> <C 1 and is completely localized if 
g'(^) » 1 in the limit i — > oo. If we just perform the thermal average (. . .)^°^ by the Wiener 
measure of a free polymer, e~'^o then we obtain 
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- 1 



U 



(0) 



(2.9) 



where e(°) = {2- d±)/2 and 



[n/T)a^ is the short distance cutoff along the length 



of the polymer. Thus the free polymer phase is obtained for > 2 (i.e., e'^^-' < 0) if 
U < U^^ oc \e'^^'>\U'^\ For U > UjP\ the polymer is free only for 



t < ll'^ = ar 



(0) 



(2.10) 



is a crossover length above which the polymer becomes localized to the pinning potential. 



From it, we can define a localization length 



(0) 



K " 



1/2 



U 



(0) 



21.(0) 



d\ >2 



(2.11) 



which characterize the typical transverse excursion of the polymer. On the other hand, for 
< 2 where e^^-* > 0, g^'^ ^ 1 for arbitrary small U in the limit of long polymer t — > oo. 
Defining the crossover and localization lengths in the same way, we have 



(2.12) 



with the polymer delocalized only in the limit f/ — > 0. Qualitatively similar behavior is 



obtained in d i =2, with 



(2.13) 



Eqs. ( p.ll| ) - ( |2.13| ) are valid as long as l^f^ ^ a. Throughout this paper, we shall only 
be interested in this critical regime, where the scaling behavior of localization length is 
insensitive to the detail shape of Up. 

Let us now return to the problem of a randomly pinned directed polymer, i.e., with 
?7(x, z) 7^ 0, and estimate the competition between the pinning and the random potential that 
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occurs in this case. Just as the locahzation of a polymer costs entropy in the pure problem, 
here, localization prevents the polymer from seeking out favorable regions of the random 
potential r) far from the origin, and therefore leads to a loss in the random component of 
the energy that could otherwise be gained even as T ^ 0. If a polymer is localized within 
a distance 1_l about the origin, then it is consisted of a number of uncorrelated segments of 
length Z|| [l^/By/^, each of which having a free energy of order Al^^ higher than that of 
the delocalized, optimal path (see Fig. 3). Thus the free energy cost of localization is of the 
order A{1^/ B)'^^"^^/^ per length. This plays the role of the entropy cost in the localization 
of a free polymer. For a very strong pinning potential, the polymer is again always localized 
completely within the pinning well because the energy per length gained, [/, always exceeds 
the random energy cost of localization, which is of the order Uq — A{a/ B)'^^'^^/^ per length. 

The effect of a weak pinning potential can be estimated perturbatively from (Hi) as just 
described for the pure problem, except that it now depends on the realization of the random 
potential r]. The natural quantity to examine is the average energy gained by the polymer 
in the presence of the pinning potential, 

5F=(Hl), (2.14) 

given that one end is fixed at the origin (as shown in Fig. 3). As in the case of the pure 
problem, this is just the average of the accumulated return probability of the polymer to the 
origin. Since the rms displacement is Bt^ in the absence of the pinning potential C/p, then 
5Fq oc U{a/ Bt'^YH to leading order in U . To determine the effect of the pinning potential, 
it is necessary to compare 5Fq to AFq ^ At^ , the intrinsic variations in the free energy [see 
Fig. 1(b)]. The dimensionless ratio that characterizes the strength of pinning is now 
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to leading order in U/Uq, with a|| = [a/ By/'^ being the short distance cutoff along the length 
of the polymer, and 

e{d) = l-d^C-d = '^-{d+l)C{d). (2.16) 

The critical dimension is (i = 1 + 1 since C(2) = 2/3 exactly 0,|l^. If <^ 1, which is the 
case when e < 0, then the weak pinning potential is irrelevant in the thermodynamic limit. 
If e > 0, then qq ^ 1 and even a small U will completely change the energy landscape of 
the randomly pinned directed polymer. For example, the free energy of the pinned state in 
Fig. 3 now becomes much lower than that of the optimal path of the unperturbed system. 
The apparent difference between the problem with point disorder and the pure problem 
is the behavior of qq at the critical dimension where e = 0: For the pure problem, (79°^ 



diverges logarithmically according to Eq. ( |2.9| ), indicating the marginal relevance of the 
pinning potential. But with point disorder, 

^0 oc U/Uo oc Ua/{AB) (2.17) 

remains finite. This naively suggests the irrelevance of a small pinning potential at the 
critical dimension. If this is true, then there will have to be a depinning transition at 
the critical dimension d = 1 + 1, since a strong enough pinning potential still produces a 
localized state. The simplest phase diagram in this case [Fig. 4(a)] will look quite different 
from that of the pure problem (Fig. 2). Fig. 4(a) is actually consistent with early numerical 
results of Kardar and Zapatocky and Halpin-Healy [^, who find a depinning transition 



in (i = 1 + 1. However, through a systematic renormalizat ion-group analysis presented 
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in Sec. IV, we will find that U itself is renormalized and diverges logarithmically, while 
all other parameters such as A and B only suffer finite renormalization. These results 
lead to a logarithmically diverging g and hence the marginal relevance of pinning at the 
critical dimension. Thus the phase diagram for the randomly pinned directed polymer 
is like that of the pure problem (Fig. 2), but with a shifted critical dimension and e, as 
shown in Fig. 4(b). This result is supported by the recent large scale numerical studies 
of Refs. |2^j3^. Renormalization group analysis can then be used to obtain the critical 



behavior at the depinning transition, in particular, the divergence of the correlation length, 

l± = a{ — -l) , (2.18) 



where i'± is the liberation exponent. 
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III. ANALYSIS IN REPLICA SPACE 



A convenient way to treat one dimensional objects like the directed polymer is the transfer 
matrix approach described in Sec. II. The full Boltzmann weight (or the restricted partition 
function) W{'x, t) of a polymer propagating from (0, 0) to (x, t) in a random medium //(x, t) 
is described by an equation analogous to Eq. (|2.8| ), except with f/p(x) replaced by f/p(x) + 
?7(x, t). All physical properties follow from disorder averages of the free energy F(x, t) = 
—Tin l^(x, t) and derivatives thereof. One approach to performing this average is the replica 
method, which exploits the identity In (x, t) = lim„^o ^(W^"(X) t)^^)- Up to an exchange 
of the n —>■ and the thermodynamic limit (see below), F(x, t) is given by the n*^ moment 
of W, Vr"(x,t) = W(x, ...x; t), where 

>V(xi,...,x„;t) = n / P[|Je-^"/^, (3.1) 

q; G {1, . . . , n} is the replica index, and Tin denotes the replica Hamiltonian 

Hn = E//^|^ (^) +UAU^)]-E p'^lUz) - . (3.2) 

For analytic simplicity, we approximate the columnar pinning potential by a 5-function, i.e., 
t/p(x) = —u6'^^{'x) with u ~ Ua'^^, from here on. Applying the transfer matrix approach 
directly to Eq. ( ^.Ij ) leads to the Schrodinger-like equation 

r^>V(xi, X,; t) = -n{n) >V(xi, x„; t), (3.3) 

where 

" f A 1 

nn) = E 'Y^l - ^^'H^.) - E r^'^i^^ - X/.) > (3-4) 

a=l { p^a ) 

As usual, Eq. ( p.3|) can be solved with the Ansatz 
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>V(xi, x„; t) = Y^ a,e-^^(")*/^<l>,(xi, x„) (3.5) 

j=0 

and 'H{n)^j = Ej{n)^j. In d± = 1 and for Up = 0, the ground state wave function is easily 
found by using the Bethe-Ansatz ^o{xi, x„) oc exp{ — J2a<i3 — xp]}. A similar 



solution follows for the ground state of a randomly pinned directed polymer confined to the 



semi-infinite plane x > and subject to an attractive potential f/p at x = [|1T|. Upon 
changing the strength of the attractive potential, one obtains a depinning transition for 
the half-plane problem when the ground state energy becomes larger than that for Up = 0. 
However, searches for a similar solution to the Hamiltonian ( |3.4| ) have not been successful. 
A naive application of the perturbation theory for small Up starting from the wave function 
^o{xi, fails likewise, because the low energy excited states of H{n) are not uniquely 



defined in the limit the number of particles n |12]. This ambiguity, arising from the 



exchange of the — > and the thermodynamic limits, is a well-known "trouble spot" for 



the replica method, and often calls for elaborate schemes with broken replica symmetry [|13 
For the problem at hand, the ground state itself does not involve replica symmetry break- 
ing However, it is demonstrated in Ref. [|12[ that states with broken replica symmetry 



could have energies arbitrary close to the ground state energy Eo{n). These will dominate 
in any perturbative calculations. 

We shall circumvent the problem of replica symmetry breaking by introducing a com- 
pletely different method in Sec. IV. For now, we consider another limit, A ^ 0, where 
the system almost decouples into n one-particle problems. Here an attractive potential Up 
always gives a bounded ground state in d± = 1. In term of the reduced length r = x/xq 
where Xq = T^/(2-u/t), the reduced Hamiltonian h{n) = 'H{n)/^^ becomes 
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(3.6) 



We now attempt to treat the case of weak inter-replica interaction A ^ uT by using the 



Hartree approximation [Q. Making the Ansatz for the wave function $(ri,r2, . . . ,r„) = 
4>i{ri)(f)2if2)---<Pnifn) and minimizing the reduced energy S{n) = J dri . . . drn^h{n)^ under 
the conditions / dr(j)'^{r) = 1 which we impose by Lagrangian multiphers Sa, we get a set of 



Hartree equations |Q . These simphfy to a single equation if we impose symmetry between 
the replicas (pair) = (j){r) and 6^ = s for all a's, 



5^ 



+ 25(r)+4(n-l)^7o 



1 J,2, 



r) — e 



0, 



(3.7) 



where 



^0-^ = A/(nT). 



(3.8) 



This equation can be solved easily with the boundary condition (j){r) = (f)'{r) = for |r| oo, 



giving a localized wave function 



^+,y^_2(n-l)(7oV'(r) 



with a negative ground state energy 



n 



where e = [1 + {n — l)gQ ^]^, and the localization length is 



-n 



l + (n-l)(7oi + i(n-l)V 



l± = xo/^/e (X (l + {n-l)gQ'^^ . 



(3.9) 



(3.10) 



(3.11) 



Note that the Hartree results Eqs. ( |3.9| ) - ( p.ll| ) are in fact the exact solutions of the 
Hamiltonian (|3.4|) in the limit 0, where the very strong pinning potential always 
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localizes the directed polymer within of the origin. However, we see that in the limit 
n — i> 0, the random potential characterized by tends to increase the localization length. If 
we extrapolate the Hartree result Eq. ( |3.11| ) to finite value of g^^, then we find an instability, 



i.e., /_L — *■ cxo as g^^ 1^. An obvious interpretation of this instability is the occurrence of 
a depinning transition. 

On the other hand, an analogous Hartree calculation for the problem of an attractive 
potential in the half plane yields again the solution ( p.9|) , except with e = (l + 2(n — l)gQ^y. 
This result suggests a delocalization transition at g~^ = 1/2 with l± oc {go~gc)~^, or z/_l = 1 
for the half-plane problem. However, the exact Bethe-Ansatz calculation |ll|] gave g^^ = 1 



and i'± = 2. Thus we see that while the Hartree solution gives the qualitative effect of a 
background random potential — a tendency to depin the directed polymer from the origin 
— it overestimates the influence of randomness and cannot be used reliably to characterize 
the depinning transition quantitatively. Indeed, we will show in the following section that 
the directed polymer is always pinned to the origin for all finite values of go in (i^ = 1 
dimensions. 

[ Note: Although our Hartree equation ( p.7|) agrees with that considered previously by 



Zhang ||3^ if we set m = 0, the physics of our solution is completely different: In the limit 
n — > which we consider, the disorder leads to a repulsive interaction gQ^5{xa — xp) between 
the replicas, which makes a change of the roughness exponent C, from 1/2 for free polymers 
to the larger value 2/3 possible. It is the pinning potential Up that confines the solution close 
to the origin. In contrast, Zhang considered the large-n limit where the interaction between 
the replicas is attractive, leading to the unphysical solution of an exponentially decreasing 
wave function even if m = 0. ] 
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IV. THE RENORMALIZATION-GROUP ANALYSIS 



As we have seen in the Sections II and III, both naive scaling arguments and the un- 
controlled Hartree solution of the replicated system suggest that the polymer undergoes a 
dep inning transition in the critical dimension d = 1 + 1, consistent with the phase diagram 
sketched in Fig. 4(a). However, we have not been able to construct a contro/Zed perturbative 
study for weak pinning potential Up within the replica formalism, despite the knowledge 
of the exact ground state. The technical problem encountered is the occurrence of replica 
symmetry breaking for the low energy excited states. In this section, we shall formulate a 
perturbative renormalization group study of the depinning transition in 1 + 1 dimensions 
without ever introducing the notion of replica. We shall take advantage of the mapping of 



the directed polymer to the hydrodynamics of the noisy-Burgers' equation |[T5|-[r7|. Many 



details of the mapping have been discussed in Ref. [IT^ and some are summarized in Ap 



pendix A. Here, we will outline our approach and quote the results. A number of detailed 
calculations are given in Appendix B. 



A. Formalism 

It is useful to introduce first a formal language to characterize the properties of the 
directed polymer in the absence of the pinning potential Up. For convenience, we shall fix 
one end ^{t) of the polymer at some arbitrary point x. This is implemented by introducing 
the one-point restricted partition function, 

Zo(x,t) = J V[^] S'^m - x)e-^"/^, (4.1) 
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where the subscript is used to indicate Up = 0. All thermal averages taken with respect to 
Zo(x, t) will be denoted by (. . .)x,t. The disorder averaged probability of finding a segment 
^(to) of the polymer at position y given that the other end is fixed at (x, t) is 

gS^'^) (x - y , t - to) = {SHm-y)k,v (4.2) 

where we use the subscript t to emphasize the dependence of G^^'^^ on the total polymer 
length t. The functional form of G^^'^^ has been investigated in details elsewhere WTj and is 



summarized in Appendix A. Here, we just mention the qualitative behavior, Gf^'^\r,T) = 
{Bt'^)~^^ for r <^ Bt'^ and vanishes rapidly as r ^ Bt'^ . 

It is also useful to consider the free energy of the one-point restricted polymer, 

Fo(x,t) = -riogZo(x,t), (4.3) 



which satisfies the noisy-Burgers' equation p6 



9,Fo(x,t) = ^^V^Fo - ^{VF,f + r/(x,t). (4.4) 



The free energy correlation function is 



a(x - y) = [Fo(x,t) - Fo(y,t)]^ (4.5) 



with 

Ctir) = 2A2p for r > Bt^, (4.6) 

Ct(r) oc A2(r/5)2^/« for r < Bt^. (4.7) 

The sample-to-sample free energy variation introduced in Eq. ( p. 41 ) is just 



AFo = JCt{oo)/2 = Af. (4.8) 



Again, the scaling form for Ct is summarized in Appendix A. As we shall see, the functions 
G, C, together with a number of other distribution functions, will allow us to compute the 
effect of an additional pinning potential perturbatively. 

In terms of Zq, the full partition function Z(x, t) of the polymer in the presence of the 
pinning potential Up{^{z)) = —u6'^^{^{z)) is just 

Z(x,t) = Zo(x,t)(e-^i/^),,„ (4.9) 

where Hi = JqcIz Up as in Eq. (|2.6|) . For convenience, we again use the continuum delta 
function supplemented by a short distance cutoff a (with u = U a'^^ ) to model the pin. 



B. The Renormalization Group Analysis 

We now investigate the effect of Tii on the "bare" system Zq perturbatively. Clearly the 
analysis is complicated by the fact that the "bare" system itself is glassy and thus highly 
nontrivial. However, to understand whether the phase diagram belongs to that of Fig. 4(a) 
or Fig. 4(b), we only need to compute the marginal relevancy of the pinning potential in 
1 + 1 dimensions. Fortunately, the randomly pinned directed polymer in c? = 1 + 1 is one 
of the very few glassy systems for which a great deal is known. In particular we note that 
there are two pieces of exact information: (i) The disorder averaged thermal displacement 
is the same as that of the pure system, i.e., 

T 



{m-m?)%, = -t. (4.10) 



due to a statistical tilt symmetry [jT^,^ preserved by the point disorder, (ii) A fluctuation- 



dissipation theorem for the noisy-Burgers' equation PB|J5S| which states that 
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d^Fo{x,t)dyFo{y,t) 



5{x - y) 



(4.11) 



in the limit of large t (see Appendix A). Eq. ( [4.11|) implies that Ct{r) = {AK,/T)\r\. Compar- 
ison with Eq. (^4.7| ) yields 26 /( = 1, and A^/B oc An/T. Along with the exponent identity 
Eq. ( |2.5| ), we have 6 = 1/3, ( = 2/3, and dimensional analysis gives B oc (A/kT)^/'^. We see 
that our "bare" problem, the fixed point of the randomly pinned directed polymer, is a fixed 
plane spanned by the axis k and Ak/T, or alternatively by the parameters A oc kB'^ and B. 
Hence all disorder averaged functions depend only on these two parameters in the limit of 
large t. For example, the bare coupling constant of the problem is defined as Qq = SFq/AFq 
where AFn = At^l'^ and 



Scaling form for C^^'^^ (see Appendix A) yields bF^ oc —{u/B)t^/^^ hence Qq oc u/{AB) oc 
uT/A as in Eqs. (p.l7|) and ( p.8|) . In subsequent calculations, we shall use Qq = u/{kB^). 
Following standard renormalization group treatment, we now consider the renormalization of 
the parameters k, 5, and u due to the perturbation Tii. [Note that the effect of the pinning 
potential will obviously depend on the position of the fixed end ^{t). If ^{t) is sufficiently 
far away from the origin, the polymer will never feel the effect of Up. In our calculations, we 
shall fix ^{t) at the origin (as shown in Fig. 3) to evaluate the maximal effect of the pinning 
potential.] 

We begin with the renormalization of the stiffness k, since the presence of a pinning 
potential breaks the statistical tilt symmetry and hence the exact relation Eq. ( |4.1(j| ). Adding 
a term —h J^dz^d^/dz) to Ti. and using the transformation ^(z) ^{z) + h{t — z)/k and 





(4.12) 
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the property ri{^{z) + h(t — z)/k,z) = r]{^{z),z) in the statistical sense, we obtain the 
renormahzed free energy (with one end fixed at the origin) 

F%tJ] = F};{0;t) -^t-u f dz{5[i{z) + h{t - z)/k])^^^ (4. 13) 

ZiK Jo ' 

to the lowest order in u. Note that Fq does not depend on h due to the statistical symmetry. 



Remembering that (([^(t) — ^(0)]^))^.^f = —TdlF\h=o where ((. . .)) denotes thermal average 
using the full Hamiltonian Ti, and defining the renormalized stiffness constant from 



we have 



K-'=^-'~u^ f^Gi'''\-h{t-z)/^,t-z)\,=o (4.15) 
= ^-\l-C^go). (4.16) 



In Appendix B, we find to be a finite constant given by Eq. (|B2D. Thus, there is only a 
finite renormalization of k which could be absorbed in its redefinition. 

Next, we consider the renormalization of the amplitude of the mean square displacement 
of the polymer 

X^ = Jdr r^Gi^''\r,t). (4.17) 
Expansion to the lowest order in u yields 



X' = (i?tC)2= / dy 



Gi'''\y,t) + -{6m-y]nin,t 



{Bt^y + u r dyy^ f dzGf'^\Q, t-z;y, t) (4.18) 

J-oo JO 



where 
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Gf^'\x -y,,t- zi; x-y2,t- z,) = -- (^(^(^i) - 1/1)5(^(^2) - 2/2))^,* (4-19) 

is a nonlinear response function of the noisy-Burgers' equation. In Appendix A, we obtain an 
expression for G^^'^^ in term of G^^'^^ using the mode couphng approximation, the vahdity 
of which will be discussed shortly. Using the relations ( [A4| ) and ( [A21| ), we obtain from 
Appendix B 

B = Bil-CBgo), (4.20) 

where is another finite constant given in Eq. ( p6|) . As for k, the small correction to B 
can be absorbed in its definition. 

To check the above result in a different way, we calculated also corrections to AFq, or the 
relation ( [4.11|) , due to the perturbation Hi. Expanding the total free energy to the lowest 



order in u gives 



d.F{x, t)dyF{y, t) = ^6{x - y) - dzCi'''^ {x - y, 0; t) (4.21) 



where 



G^r\x - y, 0; ^,t) = Fo(x, t)(5[e(^)]),, + Fo(l/, t) (5[e(^)])., (4.22) 

is another three-point response function of the noisy-Burgers equation. Again, we use the 
mode-coupling approximation to relate G^^'^) to the elementary functions G^^'^^ and C in a 
simple way. Using the relations ( [A4D , ( [A15| ), and ( |A22D , we obtain 



d,F{x, t)dyF{y, t) = — [5{x -y)- gof{x, y)] , (4.23) 



with f{x,y) ~ |a; — ?/| ^ if x + y = 0, and a faster decay for x + y 7^ Thus, the 

perturbation Tii changes the form of the correlation function. However, it does not lead to 
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a smg'it/ar renormalization of AFq, which is obtained by integrating Eq. ([4.23|) over both the 
X and y coordinates. 

Finally, we consider the renormalization of u itself. This follows from the expansion of 
the disorder averaged free energy up to second order in u: 



5F{x,t)= F{x,t) - Fo{x,t) 



1 



= (^i).,-^(^f)S,f (4.24) 

It will be convenient to integrate over the fixed point x. Using the normalization of G^^'^\ 
J^oo dxG[^'^\x, t) = 1, we have 

/u"^ /■* f°° (2 1) 

dx5F{x,t) = —ut -\ / dzidz2 / dxG\ {x,t — zi; x,t — Z2)- (4-25) 
2 Jo J-oo 

Using Eqs. (|A4D and (|A21| ), we get the following correction to u 

u = u[l + Cugo^og{t/ai\)] (4.26) 

where Cu is a finite positive constant given by Eq. ( [B12|) , and ay ^ (a/B)^^^ is the short 
distance cutoff along the length of the directed polymer (see Appendix B). The loga- 
rithmic divergence in Eq. ( |4.26| ) indicates the breakdown of the small-w perturbation for 
(7o log(t/a||) ^ 1. Thus the pinning potential is strongly relevant beyond the localization 
length 

/|| = aye^/^", (4.27) 

and the directed polymer becomes localized to the line defect. 

Since the only quantity that suffers nontrivial renormalization is the pinning strength u, 
it is straightforward to form a renormalization-group recursion relation which allows one to 
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obtain information about the polymer in dimensions d > using the calculations performed 
at the critical dimensionality dc= 1 + 1. From Eq. ( |4.26| ) and with the rescaling t' = bt, we 
obtain the recursion relation 

b^ = e{d)g + C.g' (4.28) 



where e{d) = 2 — (^{d){d+ 1) as given in Eq. ( p.l6| ). Thus the coupling constant g{b), which 



is a dimensionless measure of the pinning strength at scale b, flows to large values from any 
non-zero initial value, leading to the pinning of the directed polymer to the line defect at 
large scales if e > 0, i.e., if d < 1 + 1. But for e< or d > 1 + 1, the pinning potential 
is only effective if its strength exceeds a critical value, gc = e/C^. Thus the phase diagram 
of the directed polymer is in fact given by Fig. 4(b), like that of the pure problem (Fig. 2), 
rather than Fig. 4(a) as suggested by naive scaling arguments and the uncontrolled Hartree 
calculation. 

C. The Depinning Transition 

We now investigate the critical behavior of the directed polymer close to the depinning 
transition. Right at the depinning point gc, the wandering exponent Cc and the energy 
exponent 9c are simply 

Cc = C{d) and 9c = 9{d) (4.29) 

to 0(e) since none of the parameters k, A and B pick up divergent renormalization. The 
divergence of the correlation lengths, as one approaches the depinning transition from the 
pinned side, can be read off from Eq. ( [4.2^ ). We obtain 
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/|l = ail - 1 j and = a (^^ - Ij (4.30) 

with the hberation exponents 

i/|| = l/|e|, and zy^ = C,z/|| = C(f/)/|e|, (4.31) 



again vahd to 0(e). In d = 2 + 1 where C ~ 5/8 and e = 1/2, Eq. ([4.31|) gives z/_l ^ 1.25 



which is comparable to the results of numerical simulations: z/^ = 1.3 ± 0.6 by Balents and 
Kardar ^9], and = 1.8 ± 0.6 by Tang and Lyuksyutov It will be interesting to see 



whether Eqs. ( 4.29 ) and ( 4.31|) might be valid to all orders in e (up to some upper critical 



dimensions), as they do in the case of the pure problem without point disorders ||3TI|. This 
can be directly probed numerically in d = 1 + 1 by using correlated point disorder or by 
modifying the form of Up (see Sec. V). 



The results Eqs. ( [4.29| ) and ( [4.31| ) have been conjectured earlier by Kolomeisky and Sta- 



ley [ PU| , PT| , who combined the renormalization-group flow equations for the pure depinning 
problem with r] = and the one for only point disorder with Up = in an ad hoc way. In 
particular, Kolomeisky and Straley neglected to consider the renormalization of the stiffness 
K, B, and the amplitude of the disorder potential A due to the pinning potential Up. It is 
the lack of any divergent renormalization of these quantities, as obtained through explicit 
calculations in this study, that ensures the validity of Eq. ( |4.29| ), at least to 0(e). The lack 
of such divergent renormalization follows from the large transverse wandering of the polymer 
at the depinning transition, thus diminishing the effect of the pinning potential. (Similar 
behaviors have been found for the renormalization of the surface tension of a wetting layer 
at the wetting transition pT|.) 

Balents and Kardar ||32| obtained similar conclusions, Eqs. (|4.29| ) and ([4.31|) , by general- 
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izing the directed polymer to a D-dimensional manifold, and then computing the analogy of 
the coefficients C^, Cb and C„ in a functional renormalization group (FRG) analysis to first 
order m 5 = A — D. Since the 0(e) results ( [4.29| ) and ( [4.31[ ) do not depend on the numeri- 



cal values of the C"s, the FRG approach is effective despite the large expansion parameter 
(5 = 3) for the directed polymer. 

Finally, we comment on the validity of the mode coupling approximation used in the 
evaluation of the expressions leading to Eqs. (|4.15|) , (|4.2CI|) , ([4.23|) and ([4.26|) . Due to the 



combination of a statistical tilt symmetry and a fluctuation dissipation theorem in li = 1 + 1, 
the mode coupling approach gives the correct scaling behavior for the functions G^"^'") (see 
Refs. P6|,^ and Appendix A). This is all that is needed here since, as explained above, the 



results ( [1.29| ) and ( [4.31|) are independent of the numerical values of the coefficients C"s. On 



the other hand, the mode coupling method is known to give very good quantitative results 
even for the scaling functions [^,^,^]. Thus the coefficient C"s computed in this way 
should be quantitatively accurate. 

[ Note added: Very recently, Kolomeisky and Staley modified their renormalization- 
group analysis and obtained a power-law rather than exponential divergence of the local- 
ization length at the critical dimension 1 + 1. We disagree with their conclusions and point 



out what we believe to be the key difference between Ref. [33| and the present work: The 



analysis of Ref. is a one-loop calculation with respect to the pure directed polymer. The 
effect of point disorder is taken into account by an Ansatz which makes appropriate rescal- 
ing of the coefficients. On the other hand, the analysis presented in this section starts from 
the low-temperature fixed point of the randomly-pinned directed polymer, and takes into 
account of the extended pinning potential in a systematic treatment. The two approaches 
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should eventually lead to the same scaling behavior (see Sec. IV.D). However, the approach 
of Ref. requires a careful, consistent formulation. In particular, the choice of parameters 
[Eq. (2.4)] used in Ref. implies that not all "temperatures" renormalize in the same 
way (since Galilean invariance must be respected). This introduces some ambiguities in the 
scaling ansatz used there. ] 

D. Physical Interpretations 

As shown by the renormalization group analysis, the logarithmic divergence of the ef- 
fective pinning strength u at the critical dimension is the single most important element 
leading to the phase diagram Fig. 4(b) and the exponents ( [4.29| ) and ( |4.31| ). Discussions of 



the preceding paragraphs also illustrates the robustness of these results — the same results 
are obtained from a variety of methods, some with drastic unjustified approximations, as 
long as the scaling behaviors are properly included. We shall now provide a physical pic- 
ture of the weak localization using a phenomenological scaling theory, and argue that the 
logarithmic divergence of u should indeed be expected at the critical dimension. 

In Sec. II, we established the energy gained by the directed polymer (due to the presence 
of a weak pinning potential) to be 5Fq oc u Jq dr^Br)^^^'' ~ t^^^ in 1 + 1 dimensions. This 
corresponds to describing the optimal path of the directed polymer as a generalized random 
walk (with the wandering exponent ( rather than 1/2), and then equating the energy gained 
to the accumulated return probability of the random walk. However, such a treatment of 



the optimal path is oversimplified. As discussed in detail in Ref. , while the optimal path 
in a typical sample (or a typical region of a very large sample) is unique, there is a nonzero 
probability that there exists a different path whose free energy is within 0(1) of that of the 
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optimal path. The probabihty of finding two such paths a distance A apart is p(A) ~ A 



in 1 + 1 dimension [|T^, and the "droplet" formed (i.e., the difference between the two paths) 
typically have a length r ~ A^/^, which is of 0{/S?/'^) in 1 + 1 dimension. Thus, if one of 
the optimal path encounters the origin, then the probability of having a nearly degenerate 
optimal path also encountering the origin is p(r^/'^) ~ in 1 + 1 dimensions. In this way, 
we see that the accumulated effect of the statistically droplets existing at different scales 
leads to a logarithmic divergence, i.e., dTp(T'^^^) ~ logt. This makes the columnar pin 
marginally relevant at the critical dimension, and results in a phase diagram of the type 
depicted in Fig. 4(b) rather than that in Fig. 4(a). Note that the important ingredients 
of the above argument are only the probability of droplet formation and the shape of the 
droplets. It therefore suggests the marginal relevance of an extended defect at the critical 
dimension to be a general consequence of the droplet scaling theory. 

Technically, the droplets manifest themselves in the expression describing the renormal- 



ization of the pinning strength in Eq. ( [4.25| ), a diagrammatic representation of which is 



shown in Fig. 6. The similarity of the droplet configuration in Fig. 5 and the loop diagram 
in Fig. 6 is striking. As explained in detail in Ref. [0, the distribution of droplets is given 
by the function G^'^'^\ However, the particular distribution discussed in Ref. [|T^ has to do 
with droplets with a given width A at the same vertical coordinate, whereas the distribution 
needed here is the one with a given length r at the same transverse coordinate. Arguments 
leading to the logarithmic divergence in the preceding paragraphs assumes r ~ A^/^. This 
is validated by the more detailed calculation given in Appendix B. 



In some ways, the depinning transition here is similar to the de Almeida-Thouless |42 



line in a spin glass. The columnar pin which breaks the statistical translational symmetry 
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of the directed polymer plays a role similar to the external magnetic field which breaks the 
statistical up-down symmetry of the Ising spin glass. It should be interesting to study the 
depinning transition in the replica formalism. The encounter of replica symmetry breaking 
alluded to in Sec. Ill is no longer mysterious now — it is simply how the droplet excitations 
are manifested within the replica formalism |T^,[T^. 
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V. RELATED DEPINNING PROBLEMS 



The results of Sec. IV can be easily generalized to long range correlated point disorders, 
other forms of pinning potentials, and higher dimensional elastic manifolds. Here we shall 
describe a few interesting cases to illustrate the method. 

Following Nattermann and Medina et al [Q, we consider a random potential with 



long-range correlations of the form 



ri{x,z)ri{x',z') =2A5{z- z')\x-x'f''-^^\ (5.1) 

We shall be focus on the case d = 1 + 1, where the model (|2.1|) describes the domain wall 
of the low temperature phase of a random Ising model in two dimensions. The correlator 
{ f).l\ ) extrapolates smoothly between the case of random bond (p = 0) and random field 
(p = 1) [^,0. It is known that the roughness exponent C(p) stays at its value of 2/3 for 



short range correlated disorder as long as p < 1/4. For l>p>l/4it takes on the Flory 
value of C(p) = 3/(5 - 2p). From Eq. (^^) and e = 2 - ({d + 1), we see that e = for 
p < 1/4, and e = (1 — 4p)/(5 — 2p) < for p > 1/4. Thus for p > 1/4, a small pinning 
potential is irrelevant. As one varies u, there will be a depinning transition with exponents 
z/|j = l/|e| = (5-2p)/(4p-l) andz/^ = i/yCc = 3/(4p - 1) to 0(e). If the exponent z/y = l/|e| 
is indeed exact to all orders as in the pure case, then the above expression will be valid for 
all 1 > p > 1/4. In particular we should obtain a depinning transition with the exponents 
z/j^ = z/|| = 1 for p = 1. A numerical simulation of the depinning transition of the domain 
wall of the 2d random field Ising model should therefore be an efficient way of probing the 
"exactness" of the 0(e) result. 

We can generalize the methods described in this paper to study the pinning of a 
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directed polymer by other forms of the pinning potential Up{r). If we write Tii = 
J d^^rUp{r) J^dzS'^'-^lr — ^{z)], then the lowest order correction to the free energy is 

6Fo{t) =1^ = dz J d^^r Up{r) Gi'''\r,t - z), (5.2) 

where G'-^'^^ is the one-point distribution function defined in Eq. (|4.2| ), with the scaling 
properties described in Appendix A. Suppose Up{r) has a long tail, i.e., Up = u/r^, then 
SFo(t) ~ t^~'^''. Comparing this to AFq ~ and recalling the exponent identity Eq. (p. 51), 
we find the critical dimension of Up to be 

2 = C(4)(2 + s). (5.3) 

We expect a depinning transition at finite u above the critical dimension dc, with the liber- 
ation exponent z/|| = 1/|2 — C{d){2 + s)\. 

We may also consider a pinning potential Up{r) = —u6'^^ [r — R(-2)], where the trajectory 
of the pin Ii{z) is itself an arbitrary (but fixed) function of z. For example, we may have 
R{z) ~ z''^, describing the quenched defect trajectory of a superconductor subject to random 
collision by heavy ions. In this case, Eq. (p.2|) becomes 



SFo{t) = -u /* dz Gt{R{z),t - z). (5.4) 
Jo 

If R{z) < |t — or < C, i.e., if the transverse fluctuation of the defect trajectory is 
smaller than that of the randomly pinned directed polymer, then the 2;-dependence of the 
pinning potential is irrelevant and we have SFQ{t) ~ ut^~'^^^ as before, and the defect acts 
like a straight columnar pin. However, if R{z) > \t — z\^, as for example is the case for a 
misoriented columnar pin where R{z) oc \t — z\ [HI], then G^^'^^ is sharply cutoff and 6Fo{t) 
becomes finite for large t. In this case, a weak pinning potential is irrelevant compared to 
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the random energy gain AFq ~ t^. This type of analysis can be extended to study a large 
variety of pinning potentials, including the case where tl{z) itself is the trajectory of another 
directed polymer . 

Another interesting ramification of our considerations pertains to D-dimensional oriented 
manifolds in d embedding dimensions, i.e., z G 3?^ and $,{z) G 3?"'^^. D = 1 describes a 
directed polymer or a flux line, and D = 2 describes, for example, the domain wall of a three 
dimensional ferromagnet. The free energy scales in this case with an exponent 6 = 2(+D—2, 
which is the D-dimensional analog of the exponent identity (|2.5|) . Following Ref. p9 



we also generalize the extended defect to be ra-dimensional, i.e., Up{r) = — u5'^~"(r), and 
Hi = J d^z Up[^{z)]. The columnar pinning potential discussed corresponds to n = 1, and 
planar defects such as grain boundaries correspond to n = 2. Straightforward generalization 



of the scaling arguments of Sec. II and Eq. ( ^. 12| ) gives SFo^t) = (TYi) ~ —ut^ '^^'^ where 



t is now the linear size of the manifold. Comparing this to AFq ~ t^, we find the condition 
for the relevance of Up to be given by 

2 > C(rf)(2 + rf-n). (5.5) 

An important application here is the pinning of an interface to a planar defect, say the 
localization of a domain wall to a grain boundary in a random ferromagnet. In this case, 
D = n = 2, d = 3, and Eq. ( ^.5|) reads ({3) < 2/3. For the interface of a random bond Ising 
magnet, this is always fulfilled and the interface is always pinned by the planar defect. For 
random field systems however, ({3) =2/3 exactly in 3 dimensions, and our criterion ( p.5| ) 
is inconclusive. Straightforward generalization of the scaling argument given in Sec. IV. D 
shows that the interface is again weakly pinned by an arbitrary weak pinning potential. 



32 



ACKNOWLEDGMENTS 



We are grateful to many helpful discussions with L. Balents, D. S. Fisher, T. Halpin- 
Healy, M. Kardar, E. B. Kolomeisky, J. Krug, D. R. Nelson, and L.-H. Tang. TH is supported 
by US Department of Energy Contract No. DE-FG02-90ER40542. TN acknowledges the 
financial support of the Volkswagen Foundation and the hospitality of Harvard University 
where a part of this work was done. 



33 



APPENDIX A: USEFUL RELATIONS FOR THE DIRECTED POLYMER 

In order to make the paper more self-contained, we summarize here some useful results 
on the randomly pinned directed polymer obtained in earlier publications The 



restricted partition function of a directed polymer of length t with one end ^(t) fixed at x is 

Zo(x, t)=J [x - |(t)]e-^°/^. (Al) 

The free energy Fo(x, t) = — TlogZo(x, t) fulfills the noisy-Burgers' equation p6H3^ 

9tFo(x,t) = L^^Fo-^{VFof + r]{^,t). (A2) 

In Sec. IV, we introduced various correlation and distribution functions to evaluate the 
effect of the pinning potential Hi. To begin with, the one-point distribution function is 



miz)-y])^^, = Gi'^'\^~y,t-z). (A3) 

Again, (. . .)x,t denotes thermal average using the partition function Zo(x, t) and the subscript 
t is used in G^^'^^ to denote the explicit t-dependence of the distribution function. By simple 
scaling and normalization requirements, we have 

Gi'^'\r,T) = (Br'^r'-g./^ir/BT^), (A4) 

where gt/ri^) is finite, gt/ris) decreases sharply for s ^ 1, and / d'^-^s gt/ris) = 1 for all t/r. 
The rms displacement Xf{T) is given by the second moment, 

X^{t)= J d''^vr^G^t'^\r,T) 

= B't'^ J d'-s s%/As), (A5) 



34 



which is a weak function of t/r. It is found numerically |ll9| , |20[] that Xj^(r)/r^^ is the same 
order of magnitude for t = t and t ^ t. We fix B by the rms displacement of the free end, 
Xf{t) = BH^'^ (see Eq. This is accomphshed by choosing 



j d^^s s%{s) = l. 



(A6) 



Other than G^'^'^\ we are also interested in the higher-order distribution functions, 



T 



(A7) 



,(1,2) , Xi + X2 _ ti + t2 



*^[il,i2]*^^^ -X2,tl -t2 



y, 



Fo{x^,h){6[az)-y]).,M+Foix2.t2){6[az)-y]).„t, , (A8) 



as well as the free energy correlation function |H9 



C[ti,j2](xi - X2,tl - t2) = [Fo(xi,ti) - Fo(x2,^2)]^ 



(A9) 



where the subscript [ti, t2] denotes the larger of ti and t2- In analogy to G^^'^\ the correlation 
function has a similar scaling form. 



Ct{r,r) = AH'%/^{r/BT'^). 

The scaling function c gives the following scaling properties for C, 

= 2AH'^'^ Bt^ < Bt^ < r, 

oc 2A'^{r/Bf^/< Bt^ <^r<^ Bt^, 

For convenience, we also write the correlation function as 



(AlO) 



a(r,r)<! 



(All) 



C[ti,t2](xi — X2, tl — ^2) 



gIJ'^^O, 0) + Gg'')(0, 0) - 2G[J5)](xi - X2, ti - t2) 



(A12) 



35 



where 



G[°'J^](xi - X2,ti - h) = Fo(xi,ti)Fo(x2,t2). (A13) 



The scahng properties of Gf'''^\r,T) is easily obtained from Eq. ( JAIID . For example, from 
Eq. (^) and Cf'^V ^ oo,r) ^ 0, we have Gi°'^\o,0) = ^^(oo) = AH^'^. Below, we 



shall provide the approximate forms for all of the G'^'^'^^'s. 

It was shown in Ref. [|l^] that G^"^'"'^ can be obtained simply by adding a source term 



J(x, t) to right hand side of the equation of motion (^) and then taking appropriate 
derivatives, i.e.. 



G^-'^^ = ^ . . . -j-^ r[Fo(xi,ti) . . .Fo(x.,t„)]. (A14) 

dJ{yi,Zi) dJ{ym,Zm) 

In the context of the stochastic hydrodynamics of the noisy-Burgers' equation, the above is 



nothing but the generalized response function. In Ref. it was shown that the nonlinear 
response function G^^'^) gives the statistics of the rare but singular "droplet excitations", 
which are connected to replica symmetry breaking in the replica formalism P^JT^ . Here 



we encounter them again in perturbation theory, as we already did when using the replica 
formalism in Sec. III. 

However, unlike Sec. HI where we failed to develop a perturbative expansion due to 
the lack of knowledge of replica-symmetry broken excited states, here we can construct the 
forms of the nonlinear response functions G*^™'*^^ rather straightforwardly by exploiting a 
Fluctuation Dissipation Theorem (FDT), which the equation of motion ( |A2| ) satisfies in 
1 + 1 dimensions P^ , |iB[| . For example, the FDT gives 



l^G^^ (X - r) = ^Gf (X - y, r), r > 0. (A15) 
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Taking the limit r ^ in Eq. ( |A15|) and using the definition of G^^'^\ we immediately 
obtain 



An 



d^Fo{x, t)dyFo{y, t) = —6{x - y). 



(A16) 



Similarly, by integrating Eq. ( [A15| ) and using the scaling forms for G^^'^\ we can recover the 
scaling properties of G given in Eq. (All). 



The combination of the FDT and a Galilean invariance (corresponding to the statistical 



rotational symmetry) in 1 + 1 dimensions allows one to use a mode-coupling scheme M^PS 
to obtain the forms of the functions G^™'"). In particular, G^^'^^ and G^^'^^ are given, to a 
very good approximation, by the following set of self-consistent integral equations. 



Gf^'\x-y,t-z) =G'^'''\x-y,t-z) 



dx' dy 

^[°i^t2](^l ~ ^2, tl - t2)= '5[°'^j2](Xi - X2, ti - t2) 

\ /"OO 



1 

I2 / 

d d „{o,2)f / / ,/ d ^(1,1). , , X 
T^l^^t' -y^i -^)jr-,G,. {y-y,z-z) 



+— / dx'dy' / dt' / rfz'G(^'^)(a; - x' ,t - t')—G[r'(x' - y\t' - z 

J~oo Jo Jo ox' 

(A17) 



dy' 



dx'idx'2 / dt'i j dt'2G\^' ( 
Jo Jo 



Gt2^\x2 — X2, t2 — t'2 



dx',dx'2 ^'^'''^"^^^ ' 



where 



exp 



K T 



r > 



is the "bare" response function, and 



(A18) 



(A19) 



dx' / dt'G^^'^^ (xi - x', ti - t')G^^'^^ {x2 -x',t2- t') 

-00 Jo 



(A20) 
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is the "bare" correlation function. The mode-couphng equations (|A17| ) and ( |A18|) can be 
solved by using the scaling forms ( |A4| ) and ( |A11| ) for G^^'^^ and G^^''^\ The scaling func- 
tions g and c obtained in this way are in very good agreement with those from numerical 
simulations [|T^^^. It is found that ga{s) is approximately a Gaussian with a weak a 
dependence, and the width of the "Gaussian" is fixed by the condition ( [A 61 ) to be 1. 

The functions G^^'^) and G^^'^) can now be used to construct higher order response 
functions via the mode-coupling scheme. For instance, 



Gf''^\x-yi,t- zi, x-y2,t-Z2) 



1 



K 



/ dx' / dt'Gt^'^^ (x — x' ,t — t') 

J-oo Jo 



—G\, '{x -yi,t -zi)—G\, '{x -y2,t - z^) 



dx' 



(A21) 



'^[tiM (^^i ~ ^2, ti-t2, — y, — z) 



/ dx 



\ [ dt[Gi^'^^ {xi - x[,ti - t[) 
Jo 



dx'i 



+permutationof [(a;i, ti) ^ (3^25^2)]- 



(A22) 



The validity of the above mode coupling approximation for G^"^'^^ was discussed in detail 
in Ref. |T^. Eqs. (|A21|) and (|A22|) should capture the key scaling properties but may not 
be quantitatively accurate. We shall nevertheless use the above expressions to evaluate 
the integrals obtained in the RG analysis of Sec. IV. As our main concern is the existence 
of logarithmic divergence in the renormalization of various parameters, rather than the 
numerical value of any particular integrals, the use of the mode coupling approximation 



should be adequate. 
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APPENDIX B: RESULTS OF CALCULATIONS 



In this appendix, we compute the perturbative effect of Hi on the parameters k, B, 
and u, by using the expression for ^'^'"'"^'s obtained from the mode couphng approximation 
(see Appendix A). We are particularly interested in the t-dependence of the renormalized 
parameters k{t), B{t) and u{t) in the limit t —>■ oo. 

We start with the renormalization of the stiffness k. From Eq. (f4.16 ) and the scaling 
form (|A4|) for G^^'^\ we have 

—1 ~i f' dz it - zf ^ 



1 u 



dag'l^iO), (Bl) 



K K^B^ Jo 

where 'primes' indicate derivatives of g and = 2/3 in = 1 + 1. Thus we obtain Eq. ( f4.15| ) 
with go = u/{kB^) and 

C^ = f^dag'i,^{Q) (B2) 

which is finite. 

Next we consider the renormalization of the transverse wandering coefficient B. From 
Eq. ( [4.181 ), we have 



with 



Ib= u [ dyy^ [ dz G[^'^^ {0,t- z]y, t) 

J-oo Jo 

= — — dz dt' / dx' / dy y^Gf^'^\x' ,t — t') 

K Jo Jz J ~oo J —oo 



^G^^^'\x'-y^t')^G^^^'\x\t'-z)^ (B4) 
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where we used the mode-couphng approximation Eq. (|A21|) for G^^'^\ Using the scahng 



form (|A4|) for G^^'^\ and noting that g{s) is symmetric in s, we find 

Ib = -iBH'"^)2goCB (B5) 

where 

Cb= ^ da dsi-s)g[,_^yAs) ~9i,r{sa^ /r^). (B6) 

JO JO J — oo 



Again is finite since g are normahzed and sharply cutoff for large argument. We thus 
obtain the result Eq. ([4.20| ). 



Finally, we consider the renormalization of u. The expression given by Eq. ( [4.25| ) can 



be described diagrammatically as in Fig. 6. If we use the mode-coupling approximation 
Eq. ( [A21| ) for G*^^'^-', and note the normalization condition / dxG^t'^\x,T) = 1, we find 

6F{t)=—ut + — / dzidz2 I dxGf^'^\x,t — zi, x,t — Z2) 
2 Jo J-00 

= -ut-^ fdt' f dz, f dz2 r dx'^G^l''\x',t' - z,)^G^l''\x',t' - Z2). (B7) 

ZK jo jo jo J — cxd ux Ox 

Using the scaling form for G'^^'^^ again, we obtain 

SF{t) = -ut--— / dt'hit'), (B8) 
ZkB'^ Jo 



where 



4(r) = r r ds t,tM K/tM/i)- m 

JO ti 12 J-00 



We shall see that /„ is actually divergent. To regularize the integral, we insert a ultra-violet 
cutoff scale ay oc (a/i?)^/^, since in our model, the columnar pin Up is really a potential well 
of finite size a; it is only approximated by a delta function at scales much larger than a. 
Eq. ( P9|) then becomes 
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4(r)= / ^F(ti/r,a||/r), (BIO) 

Jay ti 

where 

•J Q, I L -J — OO 

Clearly the divergent part of /u(t) comes from the limit T/ay — > in Eq. ( piO|) . So to 
leading order, we have luij) = log(T/a||)C„, where 

Cu= lim -F(a, a) 

/CXD ^/-r POO 
^ J_J'9'oo{s)9'oo{s/cr'^), (B12) 

which is positive definite since = 4/3 > 1 and g is well behaved. The effective pinning 
potential u can now be defined as 

u = -^^SF{t) = u[l + goCJogit/aii)], (B13) 

which is quoted in Eq. ( |4.26| ). 
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FIGURES 

FIG. 1. (a) The solid, dashed, and dotted paths are the optimal paths of a randomly pinned 
directed polymer with one end fixed at xi, X2, and X3 respectively. The rms displacement of 
the free end is X{t) oc t^, where t is the polymer length. For small separation of the fixed end, 
X2 — xi <C X(t), the two optimal paths typically merge at a length r cx |x2 — a^il^''^ from the 
fixed end. For much larger separations, X3 — xi 3> X{t), the optimal paths become completely 
uncorrected, (b) Variations in the directed polymer's free energy as a function of the fixed end 
position, X. For large separation of the fixed end, typical free energy differences are AFq oc t^. 

FIG. 2. The phase diagram of a pure directed polymer pinned by a short-ranged columnar 
pinning well Up. The arrows indicate the effective pinning strength U in units of the entropy cost of 
localization, Uq'^ = T^/(2Ka^). The directed polymer is pinned in d± < 2. In > 2, a depinning 
transition at U = U^^ oc {d_\_ — 2)U^^ (indicated by the asterisk) separates the pinned and depinned 
phases. 

FIG. 3. The solid line depicts a directed polymer localized to a columnar pin at the origin. 
The localization distance is The polymer is consisted of a number of uncorrelated segments of 
length Z|| oc . An optimal path far away from the origin is depicted by the dashed line. 

FIG. 4. (a) A possible phase diagram for the directed polymer pinned by both the columnar 
pin and a random background. The arrows indicate the effective strength U of the columnar pin, 
in units of the random energy loss, Uq = A{a/ B)^^~^^/^ . The critical dimension is d_L = 1. Unlike 
Fig. 2, this phase diagram contains a depinning transition in all dimensions, (b) An alternative 
phase diagram having the same structure as that of the pure one (Fig. 2), but with a shift in the 
critical dimension to = 1. 
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FIG. 5. The solid line is an optimal path which intersects the origin at z = zi. The dashed line 
is a path that is nearly degenerate in free energy. It intersects the origin at z = Z2- The difference 
between the two paths is called a "droplet". It has a length r and a width A ~ r^. 

FIG. 6. A loop diagram describing the renormalization of the pinning strength u given in 
Eq. (4.25). The three-point function G^^'^^ is given by the convolution of a vertex (the black triangle) 
and three G'^^'^^'s ( lines with arrows). The mode-coupling approach Eq. (A21) approximates the 
vertex by the "bare" vertex, i.e., two spatial derivatives (see Appendix A). 
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